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Abstract 

In this paper we show, in a systematic way, how to relate the Kepler 
problem to the isotropic harmonic oscillator. Unlike previous approaches, 
our constructions are carried over in the Lagrangian formalism dealing, 
with second order vector fields. We therefore provide a tangent bundle 
version of the Kustaanheimo-Stiefel map. 



1 Introduction 

Reduction procedures, providing a way to link a given dynamical system with 
one on a lower dimensional manifold, have been extensively studied as an help 
in integrating the dynamics (pQ, 0, 0, EQ, It has been noted, however, 
that often a sort of inverse of reduction, an "unfolding" procedure, can be more 
effective for this purpose ( 7 ). It has been shown (J5j) that different classes of 
completely integrable systems arise as reduction of free or "simple" ones with 
higher degrees of freedom: in this approach, one can obtain the solutions of the 
given dynamical system as a projection of the ones of the higher dimensional 
system (exactly known). The projection method would be therefore of interest 
in the problem of integrating the dynamics to develope an unfolding tecnique. 
It is clear that there is no unambiguous way to recognize that a given system 
can be obtained as reduction of a free one, since there are arbitrary elements in 
the reduction procedure itself (essentially the choices of the submanifold of the 
carrier space invariant for the dynamical system and of the equivalence relation 
on it, cfr. [S], [3]). Taking this into account, one can try to focus on the 
main aspects of a dynamical system, so to impose some conditions that the 
related system has to satisfy and in this way reduce ambiguity. In this paper 
we develope such an approach for the unfolding of the Kepler problem in three 
degrees of freedom. This system has been widely studied ( |TJ]], QU, [*H3). 
|14|. |15|) and its relation with the isotropic harmonic oscillator in four dimen- 
sions through the Kustaanheimo-Stiefel map is well known ^B]. We recover this 
relation in a more general way, this may be possibly used for the unfolding of a 
general system. 
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Here we only consider the Kepler problem in three degrees of freedom, but it 
seems to us that our procedure can be fruitful also in higher dimensions (for a 
different approach see J7|, [TBI)- 
The paper is organized as follows. 

In section 2 we briefly recall the main elements of a reduction procedure, giv- 
ing two examples of completely integrable system arising as reduction of a free 
system. 

In section 3, after recalling few elements of Lagrangian formalism, we point out 
the main aspects of the Kepler problem in three dimensions, since they will be 
relevant in our unfolding procedure; moreover, we introduce some properties of 
reparametrized vector fields. 

In section 4 we recover, from general consideration, the Kustaanhcimo-Stiefel 
map as the main tool in our unfolding procedure. 

In section 5 we arrive at an unfolding system for the Kepler problem, and rec- 
ognize that it coincides, up to reparametrization, with a family of harmonic 
oscillators. 

In section 6 we characterize the symmetry of the unfolding system and find out 
the subalgebra of the constants of the motion for the Kepler system. 



2 Reduction procedure 

In this section we point out the main characteristics of a reduction procedure for 
second order systems, since a thorough understanding of reduction procedure is 
essential in performing its "converse" , the unfolding procedure. In particular, 
we treat two cases of nonlinear (completely integrable) systems arising as re- 
duction of free ones in higher dimensions, in order to provide some motivations, 
through explicit examples, for an unfolding procedure. 

Let r be a dynamical system on a carrier space M., a smooth manifold, namely 
a vector field in x(-M) ( we don't consider additional structures). Any reduction 
procedures requires two ingredients: 

i) a submanifold E of M. which is invariant for T, that is to which T is 
tangent: 

r(m) S T m E VmeS (2.1) 

ii) an equivalence relation wonS which is compatible with T, that is (calling 
$^ the flow of r): 

ra«m'^$[.(m)Ri${.(m') Vm,ra'eE (2.2) 

The reduced carrier space is the quotient manifold £ = £/ rs and the reduced 
dynamical system T is the projection of Ts (i.e. T restricted to E) on E with 
respect to the natural projection 7r:E^E/«h m — > m — [m]~ (well defined 
thanks to the compatibility condition). 
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We recall that a vector field f € %(£) is said to be the projection of € x(£) 
with respect to tt, or that is projectable onto T with respect to tt if the 
following diagram is commutative: 

TE^-TS (2-3) 



rv 



£ 



Following this line, we give two examples of reduction of free system that gives 
rise to nonlinear completely integrable systems; for further details see 0], [5], 

0, El 

2.1 Example: radial reduction 

Let us consider a free system in 3 dimensions, described by the following vector 
field on TR 3 : 

r = «^ i = 1...3 (2.4) 
corresponding to the following equation of the motion: 

f = (2.5) 

where r is a vector in R 3 . 

We will obtain a submanifold £ in TR 3 fixing the value of the energy, £ is 
invariant since E is a costant of the motion; the equivalence relation on £ will 
be given in terms of the group 50(3) of symmetry for the system. 
Introducing polar coordinates r = rr , one has 

r = fr + rr ; r = rr + 2ff + rr (2-6) 

and, using the vectorial identities: 

f-f=l; f-r = 0; (r) 2 = -f ■ r (2.7) 

and the equations of the motion, one gets: 

f = -rfr = r(f) 2 (2.8) 

From the first eauation l2.tjl it follows that: 



T r 

so Eq. 12.81 becomes: 



(r) 2 (r) 2 

f = ^L-^L (2.io) 



r r 
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For a free particle 2E = (r) 2 ; on the submanifolds of fixed energies 

E £ = {(r , v) e TR 3 h (v) 2 = 2E = cost} (2.11) 

one obtains: 

r = ^- (2.12) 

r r 

Since the equation defining is invariant under the action of 5*0(3), the action 
of this group on TM. 3 can be restricted to an action on E#. This provides an 
equivalence relation « SO ( 3) on Ej: two points of E# are equivalent if they are 
connected by a transformation of SO(3). So we obtain the quotient manifold: 

E £ = E £ /S0(3) (2.13) 

The vector field T restricted to E# is projectable since Eq. I|2.12|l is rotationally 
invariant; one has: 

~ d (IE (ff\ d , n , s 

T = v r —+( — U- 2.14 

or \ r r J ov r 

Remark: if one considers the submanifolds of fixed value for the angular 
momentum £, from £ 2 — (r) 2 (r) 2 — (r • r) one obtains: 

f -v — — — (2 1 

r dr r 3 dv r 

2.2 Example: the Calogero-Moser system 

We consider again a reduction of the free motion in R 3 ; if we parametrize the 
elements of M. 3 using 2x2 symmetric matrices: 

the equations of motion for the free particle in 3 dimensions can be written as: 

X = (2.17) 
and the Lagrangian function becomes: 

Lo = \^X 2 = \{±l+±l+±i) (2.18) 

The matrix: 

M=[X,X] (2.19) 
is a constant of the motion; being antisymmetric, it can be written as: 

M =£a, CJ = ( \ o ) ( 2 ' 2 °) 
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where £ is the modulus of the angular momentum. 

Since X is a real symmetric matrix, it can be diagonalized by elements of the 
rotation group SO(2): 

X = GQG^ 1 (2.21) 

where: 

Q=( q Z ° ), G=( C ° SIP 8[nip ) (2.22) 

\ q2 J V — sin <f cos <p J 

From this one obtains: 

X = [GG- 1 , X] + GQG- 1 = G ([G" 1 G , Q] + Q) G" 1 (2.23) 

and: 

M= [X,[GG-\X\] =G[Q, [G _1 G, g]]G _1 (2.24) 

Moreover: 

G- 1 G = G G- 1 = <pa (2.25) 
Using these relations one can evaluate the trace of Ma to find the value of t: 

e=-±Tr(Ma) = <p(q 2 - qi f (2.26) 
from which one can obtain <p. Deriving 12 . 231 with respect to time, one gets: 

X = ip 2 [a, [a,Q]] (2.27) 
that is, substituting (p from Eq. 12.261 

2£ 2 



qi 



(?2 - qi ) 3 



2£ 2 

h = 7 ^ (2.28) 

(92 - qiY 

Since I is a constant of the motion, the condition £ = const provides invariant 
submanifolds T,£ in the tangent bundle of the space of symmetric matrices; since 
this space covers TR 2 , one can project from each one of these submanifolds and 
find a family of dynamical systems in the variables 92, that is a second order 
dynamical system on TR 2 . When £ = we get the free particle in 2 dimensions; 
when I ^ we get the Calogero-Moser system in R 2 ([F^, [H],|H]). 



It is clear from these examples that having reduced the dimensions of the dy- 
namical system does not always provide a simpler system: we have obtained 
nonlinear systems starting with a linear one; moreover, in the first case the 
integration is not immediate. Such cases suggest that it could be reasonable 
for the integration of the dynamics to investigate if a given nonlinear system 
could arise as reduction of a linear one, so to obtain the solutions of the former 
projecting those of the latter. 

We shall try now to apply an unfolding procedure to the three-dimensional 
Kepler problem. 
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3 The Kepler problem 



Our first goal is to enumerate the main aspects of the Kepler system, in order 
to understand what properties one requires for the unfolding system. 
Before doing this we recall some elements of the geometry of Lagrangian for- 
malism on the the tangent bundle, assumed to be the carrier space for the 
dynamics. 

3.1 Few elements of Lagrangian formalism 

It is well known that it is possible to formulate the Lagrangian formalism in an 
intrinsic, coordinate free version, involving objects characteristic of the geom- 
etry of the tangent bundle; this is what we briefly recall here, referring to |21j 
for details. 

On a tangent bundle there are two natural tensor fields, they essentially char- 
acterize its structure: the vertical endomorphism S and the dilation vector field 
A, given in natural coordinates by: 

S = dx ^h A = ^ (3 - 1} 

In terms of these objects, a second order vector field T can be defined in intrinscc 
terms by: 

S(T) = A (3.2) 

Moreover, it is possible to associate to S a generalized derivation d$f = d/o 5 1 , 

in coordinates dsf — -—^dx l . 

ov l t 

Given a Lagrangian function £, one can define the Cartan 1-form of &: 

9c := d s H (3.3) 

and the Cartan 2- form: 

lol = ~d8 L (3.4) 

If loz is non degenerate, i.e. symplectic, L is said to be regular. 
Having introduced these objects, it is possible to write the Euler-Lagrange equa- 
tion in the following coordinate-free way: 

C T e L - d£ = (3.5) 

where T is a vector field to be determined. When H is regular, T is second 
order: in this case, after introducing the energy function 

£ L :=C A L-L (3.6) 

one can rewrite equation 13. 51 in the following equivalent way: 

i r uJc = d£ L (3.7) 
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So, when the Lagrangian function is regular, we have obtained a symplectic 
formulation of the Lagrangian dynamics directly on the tangent bundle. It is 
now possible to define Poisson brackets on the tangent bundle as one usually 
does on any symplectic manifold: 

:=w(X /t X fl ) (3.8) 

where the vector fields are obtained solving ix,W£ — df and iy^ g ioa = dg. 

In this way the vector space of functions on the tangent bundle acquires the 

structure of a Lie algebra under Poisson bracket. 

We note that now the symplectic 2-form depends on L and so does the Poisson 
bracket (see [2*2*]). 



3.2 The Kepler system 

By Kepler system we mean the model of a point particle subjected to a central 
force proportional to the inverse of the square of the radius: 

F* = -£ (3.9) 

The configuration space of the Kepler system is R 3 minus the origin; in the 
following we will denote R™ — {0} = Rq for brevity. 

To the equations of the motion, that is a second order system of differential 
equations in R 3 , it is possible to associate a first order system on the tangent 
bundle of R 3 , whose solutions are the integral curves of the vector field 1 : 

_ i_9_ kx^d_ 

K dx % r 3 dv l 

The vector field IV € x(TR 3 ) is a second order vector field (cfr. sec. 13. If) : 
it defines a system of second order differential equations on the configuration 
space. 

It is well known that the Kepler problem admits a Lagrangian formulation, with 
Lagrangian function: 

H K = \v i v i + - (3.11) 
2 r 

where r — V x i x l and we have set the mass m = 1. 

The energy associated to this Lagrangian is (cfr. Eq. (|3.6(l ): 

£^ = I V V-- (3.12) 
2 r 

and the Cartan 1- and 2- forms are, respectively: 

9 LK =v l dx l (3.13) 
= d9c K — dx* A du* (3.14) 



1 in the following latins indices run from 1 to 3, greek ones from to 3 and the summation 
over repeated indeces is adopted 
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such that one has: 

ir K UL K = d£ t (3.15) 

We notice that we can express all these objects in terms of the Euclidean metric 
ofM 3 . 

The Poisson bracket obtained from the symplectic form u>n K as in Eq. I|3.8|l is: 

df dg df dg 
= (3 ' 16) 

It is well known that constants of the motions for the Kepler problem are: 

Li = e ljk x l v k (3.17) 
At = ( — +e ijk v j L k j (3.18) 

respectively the angular momentum vector and the Runge-Lenz vector, satisfy- 
ing the equation: 

UA, = (3.19) 

We recall that in the association between symmetries and constants of the mo- 
tion by means of the Lagrangian symplectic structure uj^ only the first one arise 
from point transformations, while the second comes from a so called "extended 
symmetry" , i.e. a dynamical symmetry which does not respect the tangent bun- 
dle structure of the carrier space. 

The commutation relations among the constants of the motion are the following: 

{Li,Lj}a K — eijkLk; {Ai,Aj}c K — —22,^ijkLk 

{L i ,A j }c K = e ijk A k (3.20) 

Since the energy function £n K belongs to the center of the algebra of the con- 
stants of the motion, we can rescale the Runge-Lenz vector by appropriate 
functions of the energy (of constant sign) . We can divide TK 3 into three sets, 
the two open regions when &c K < 0, and E-a K > and their boundary £-a K = 0: 
when £.c K < 0, we can rescale A 1 by y/—28,£ K , to get an o(4) algebra; where 
&C K > 0, we can rescale A % by y / 2£,c K , to get an o(3, 1) algebra; when = 
we get the Eucledean algebra in 3D. 

The Kepler system is completely integrable and maximally superintegrable. We 
recall that a system with n degrees of freedom is said to be maximally superinte- 
grable when the largest number of constants of the motion which are essentially 
independent is 2n — 1; a set of functions f 1 ,. . . f G is said to be es- 

sentially independent if the set {x l £ M h df 1 (x) A ... A df k (x) = 0} has no 
interior points. The Kepler system has 2n-l=5 constants of the motion essen- 
tially independent, they are the components of the angular momentum and of 
the Runge-Lenz vector (only five of them are independent because of the con- 
straint enm . 

Another important feature of the Kepler system is that it has closed orbits for 
an open portion of the phase space, the one corresponding to negative energies. 
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In the search for a system the Kepler problem is a reduction of, we are sug- 
gested to impose the condition of closed orbits; moreover, it seems reasonable 
to preserve also the property of superintegrability. Along the line of section [21 
the first attempt would be in terms of linear maximally superintegrable systems 
with closed orbits, and the only system that satisfies all these conditions is the 
harmonic oscillator. 

We have required that the systems from which one can obtain the Kepler prob- 
lem by reduction have to share some properties of the Kepler problem, and in 
this way we have restricted the ambiguity in the unfolding procedure. 
Few comments are in order at this point. First, there is an obstruction to the 
correspondence between harmonic oscillator and Kepler problem given by the 
energy-period theorem. It says that, if the period T of a periodic Hamiltonian 
system is a (at least) C 1 function of the carrier space, then dH A dT = 0; so the 
period is a function of the energy: T — T(E) (see e.g. 

If we have a system T which reduces to T, as in section |21 from the equi variance 
of the flows (in the notations of section [2J: 

& t O 7T = 7T O $*, e (3.21) 

it follows that the period T of the unfolding system is also a period for the 
reduced system: 

$£+ T (m) = ¥ r (m) =*> $*+ T (m) = $~(m) (3.22) 

Thus, an isocronous system cannot reduce to an non-isocronous sytem, and vice 
versa. 

While for the Kepler problem orbits with different energies do have different 
periods, it is well known that the period of the harmonic oscillator is constant, 
independent from the energy: from the previous result, it follows that we shall 
not be able to map the isotropic harmonic oscillator motions onto the motions 
of the Kepler problem. So, if integral curves of the isotropic harmonic oscillator 
project to integral curves of the Kepler system, to match the periods we are 
obliged to use a one parameter family of different oscillators whose period de- 
pends on the energy of the Kepler system, and the correspondence can be done 
only on submanifolds of fixed energy. 

The second comment concerns another important difference between the two 
systems: while the vector field of the harmonic oscillator is complete, the one 
of the Kepler problem is not complete. So the problem arises of how to link 
two such systems. But a well known theorem (see e.g. [22], Chap.V-1.8) as- 
serts that, given a vector field X on a paracompact manifold Ai, there exists a 
strictly positive function / on Ai, of the same differentiability class as X, such 
that X = / -X is complete. Thus we should expect that at a certain point in our 
unfolding procedure there will be a need for a reparametrization of the vector 
field involved: therefore it is appropriate to point out some aspects connected 
with reparametrization. 
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3.3 Reparametrized vector fields: some properties 



First of all, we notice that the integral curves of X and X = / • X, are the same, 
but their parametrization changes. 

Anyway, the constants of the motion do not depend on the parametrization, 
since: 

Cyji = => C ± h = fC*h = (3.23) 

However symmetries, represented by vector fields, do depend on the parametriza- 
tion, since from [X, Y] = one only gets [/X, Y] = (£y/)X 
Obviously, the reparametrization of a vector field does not preserve some other 
properties. For intance, if X is Hamiltonian, from 1-^10 = dh it follows that 
ijxw = fdh, so X is not Hamiltonian unless d/ A dh = 0. 

Moreover, if X e x(TQ) is a second order vector field, X in general will not 
have the same property; in fact: 

S(X) = A =>■ S(fX) = /A (3.24) 

However, it is possible to endow TQ = M. with a different structure of tangent 
bundle with respect to which the reparametrized vector field is second order. 
We now sketch a possible procedure, referring to |2U for details. 
We have already noticed that two objects characterize the geometry of the 
tangent bundle: the vertical endomorphism S and the dilation vector field A. 
Actually, it is possible to prove that, given a vector field and a 1-1 tensor field on 
a manifold M. verifying the same joint properties of A and 5, they characterize 
uniquely the tangent bundle structure on M. , that is a unique manifold 23 exists 
such that M ~ TB (see |E2j for details). 

Using this result, given a vector field X on a manifold A4, it is possible to 
provide M. with a tangent bundle structure with respect to which X is second 
order, if some appropriate conditions are satisfied. 

Let dinxM = 2n, and let g % be n functionally independent functions on M. , that 
is dg 1 A dg 2 A ... A dg' 1 ^ 0; let us suppose that n functions /' exist such that: 

f = £xg z (3.25) 

(possibly only locally) and that they satisfy the relation: 

dg 1 A ... A dg n A d/ 1 A . . . df n ^ (3.26) 

The set (f l ,g l ) will provide a set of coordinates on M. (at least locally); so one 
can build a 1-1 tensor field and a vector field: 

S--=W®^ A: =/^ ( 3 - 27 ) 

that verify the joint properties we mentioned above. By virtue of the theorem 
in |2S], they define (at least locally) on M. a structure of tangent bundle, where 
the g % can be considered the base coordinates and the /* the fiber coordinates 2 . 

2 we notice that the g 1 are an algebra with respect to the pointwise product, while the /* 
are a module space 
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The vector field X turns out to be second order with respect to this structures, 
since: 

5(X) = A (3.28) 
using Eq. 13.251 that was our main assumption. 

If equation 13.251 is verified locally, one has to repeat this construction for each 
open set where Eq. 13.251 is verified and take into account transition functions 
which should behave as "point transformations" , i.e. tangent bundle isomor- 
phisms. 

Incidentally, we notice that the same manifold A4 can acquire different tangent 
bundle structures. For instance, one can multiply the same set of functions g % 
by a function of some constants of the motion C a for X (that is CyJ3 a = 0); 
so one has: 

G i = h{C a )g i £ x G l = h(C a )g i := F i (3.29) 

so the functions (G l , F l ) defines different 5" and A', and so a different structure 
of tangent bundle on A4 . 

4 The Kustaanheimo-Stiefel map 

So far, from general considerations, we have considerably reduced the ambiguity 
of the unfolding procedure, having identified as a candidate unfolding system 
only (a family of) harmonic oscillators. Of course, we still have ambiguity on its 
degrees of freedom: the minimal condition we can impose is that they have to 
be stricly greater then those of the Kepler problem, since the latter has to arise 
as a reduction of the former. So, calling n this number of degrees of freedom, 
we have n > 3. For simplicity one can consider first the cases with lowest 
n, and investigate if the correspondence can be realized for, say, n = 4. So 
one has to look for a map TR^J — ► TM. 3 . Since the angular momentum for the 
Kepler problem is a constant of the motion associated to point transformations, 
and since constants of the motion of the reduced sistem are projection of the 
unfolding ones, we can ask that the reduction procedure respect this structure, 
that is it comes from a map between the configuration spaces n : R„ — ► R„. 
Since we have to perform a reduction, we look for a covering of M. 3 with R*. 
Because M. 3 = S 2 x R + and R^ = S 3 x R + we may start from a covering map 
ttjj : S 3 — > S 2 and extend it. Identifying S 3 with SU(2), we can represent it in 
terms of matrices: 

S ={-y 0+i y 3 yi -iy 2 > l~ £ = detfo) = 1. (4.1) 

and so define the Hopf map ([2H1) by setting: 

JCS : s — > x h sa^s^ 1 — x l Oi (4-2) 
where cr, are the Pauli matrices (and oq = I). 

There are several ways to extend the Hopf map to R^ — > M. 3 . A natural one is 
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obtained by introducing polar coordinates in M.„ = S 3 x M + and setting 



g = Rs with s £ SU(2), ReR+ 
and, recalling that s _1 = s* for s € SU(2), we define: 



/CiS :g^a;h x k <7k = 9^9 = R 2 scr 3 s 1 



or, alternatively: 



x\ = 2(y 1 y 3 + y 2 y ) 
x 2 = 2(y 2 y 3 - yiy ) 

X 3 



vl + vl-vl-yl- 



From detg = R, we find: 



= det{ga 3 g f ) = -dct( 3 t 5 ) = -R* 



or: 



I? 2 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



From Eq. I|4.4|l it follows that the one parameter group exp(iAo3) defines the 
fibers of the fibration U(l) — > M.f t — > ; we notice that this group acts by right 
moltiplication. 

This map extends naturally to TRj -> TI 



?3. 



vi = 2(yiu 3 + y 2 u + y%Ux + y u 2 ) 
T(/C<S) : v 2 = 2(y 3 u 2 + y 2 u 3 - yiu - y ui) 
V3 = Viui + U2U 2 - y 3 u 3 - y u 



(4.8) 



and one has the fibration TU(1) -> TM^ ~> TM^. 

In this way we have recovered the Kustaanhcimo-Stiefel map, introduced in |16|. 
where the authors extablished a relation between the solutions of the isotropic 
harmonic oscillator and those of the Kepler problem. Here we have tried to 
arrive at this map in a constructive way, as far as it has been possible: we 
have made some general requirements, but also arbitrary choices, for instance 
the number of dimensions and the extension of the Hopf map; anyway, all our 
considerations relied on general properties of the system, so could be hopefully 
repeated for other systems, with the appropriate changes. 



5 Unfolding of the Kepler problem 

Having introduced all the elements we need, we can now look for the motions 
in 4 dimensions which may be related to Keplerian motions in 3 dimensions. At 
this point the Lagrangian character of the Kepler problem may be put to work. 
Since iv K ^L K — &&l k , pulling back lol k and with TICS, after appropriate 
considerations, we can obtain a symplcctic vector field whose projection is the 
Kepler vector field. 
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As for the pull back of ujc k and £-c K we may instead consider the pull back of 
the Euclidean metric dx l <g> da; 1 , from which Lk and E^ K are derived: 

(TKSYi&x 1 <g> dx J ) = (T7GS)* (^Tr[{dx k a k ) <8> (dx fc a fc )] 

= iTr[d( ff( T3g t )®d(. 9 a 3 . 9 t)] (5.1) 

Using that g = i?s, one gets: 

d(0o-35 f ) = di? 2 sazs- 1 + R 2 d(sa 3 s- 1 ) (5.2) 

so: 

^Tr(di? 2 SC73S- 1 + i? 2 d(scr 3 s- 1 )) <g> (di? 2 scras- 1 + i? 2 d(scr 3 .s- 1 )) = 

i[di? 2 <g> d J R 2 Tr(scr 3 s" 1 scr 3 s" 1 ) + i^Tr^S^S -1 ) <g> d(scr 3 s" 1 )] 
+R 2 dR 2 ® Tr[(scr 3 s- 1 )d(scr 3 s- 1 )] + i? 2 Tr[(scr 3 s- 1 )d(scr 3 s- 1 )] ® di? 2 ] (5.3) 

The trace in the first term is 1, the last two terms are zero because the trace 
involved is zero. As for the second term, using that: 

d(sa 3 s^ 1 ) = d(s)a 3 s _1 + S(T 3 d(s _1 ) = ds a 3 s^ 1 — scr 3 s~ 1 ds s _1 = 

= [s^ds, cr 3 ] = i[a k 9 k ,<j 3 } = 2ie k3l <r i e k (5.4) 

where the square bracket stands for the commutator, and we have used the 
relation s _1 ds = icr k 9 k , where 9 k are the left invariant one-forms of x*(5'^(2)). 
So we get: 

Tr^scrgs- 1 ) <g> d(scr 3 s" 1 )] = 6 1 ®6 1 + 9 2 ® 9 2 (5.5) 

As for the pull back of the metric, one finally obtains: 

(TICSy(dx l (g>dx l ) =4R 2 [dR<g>dR + R 2 (9 1 <g)9 1 +9 2 <E)e 2 )} (5.6) 

Of course, it is a degenerate quadratic form; but a comparison with the confor- 
mally flat metric: 

g c = 4i? 2 [di? ®dR + R 2 (9 1 ® 9 1 + 9 2 <8> 6 2 + 9 3 ® 9 3 )} (5.7) 

suggests that we complete it by adding the missing term 9 3 ®9 3 . We notice, in 
fact, that it doesn't affect the kinetic energy in since if one has 3 in TKf : 

T = 4i? 2 [i? 2 +i? 2 (# 2 + # 2 + fl 2 )] (5.8) 

the last term, the added one, is invariant under the right multiplication by 
exp(iAcr 3 ), the group with respect to which one has to perform the reduction; 

3 in the following 9 k = irT* 4 9 k 
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so in the reduction procedure, choosing R A 9\ = as the invariant manifold Eo 
(see section |2J) gives the Kepler system (we will show later that it actually is 
invariant). Incidentally, we note that i? 4 #3 is the Hamiltonian function for the 
tangent lift of the former U(V) action. 

With these considerations, after the pull-back we get a "conformal Kepler prob- 
lem" in M 4 , described by the Lagrangian: 

l = 2r 2 [r 2 + R 2 {ei + e 2 2 + el)] + jp (5.9) 

to which is associated the energy: 

E L = 2R 2 [R 2 + R 2 {8\ + B\ + - A (5.10) 

and the Cartan 1-form and symplectic 2-form are given by: 

9 L = AR 2 {RdR + R 2 9 k 6 k ) 
lu l = -4R 2 dRAdR- 16R 3 9 k d R A 6 k - 4 R 4 (d9 k A 9 k + 6 k d0 k ) (5.11) 

For future convenience we express all these objects in cartesian coordinates 
(y Q , u a ): 

£=4i? 2 iuV + A ; £ £=4 i? 2 i u V-A (5.12) 

9 L = 4i?Vdy i uj l = -8j/ /3 u Q dy /3 A dy a - 4R 2 du a A dy a 

The dynamical vector field obtained from it^>l — d£^ is: 



k uPy? \ f) ( 5 - 13 ) 

2i?6 y R2 U 

At this point, we have completely determined a dynamical system in R 4 that 
gives the Kepler system by reduction. More precisely, one first has to restrict 
to the submanifold So (see sec. 0) where R A 6^ = 0, which is invariant for T 
since Cr{R A d^) = 0. On this submanifold, one chooses the equivalence relation 
defined by the action of the group U(l), tangent lift of the right multiplication 
by exp(«A(73), i.e. two points are equivalent if they belong to the same orbit of 
the group. This equivalence relation is compatible with T, because the group 
under consideration is a simmetry group for T, since its Hamiltonian i? 4 #3 is a 
constant of the motion. The reduced system obtained in this way is the Kepler 
system on TM%. 

Moreover, one can consider the algebra of the constants of the motion of the 
conformal Kepler problem, where the Poisson bracket is obtained from the sym- 
plectic 2-form as previously shown: it is possible to obtain the constants of the 
motion of the Kepler problem as a subalgebra of those of F. We postpone such 
considerations to a following section. 
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5.1 Relation of the unfolding system with a family of har- 
monic oscillators 



For the reasons explained in the previous section, we look for a relationship be- 
tween this system and a harmonic oscillator, and we expect that such a relation 
involves the systems at hypersurfaces of fixed negative energies. Moreover, we 
notice that a reparametrization is required, since the vector field T is still not 
complete. 

As stated in the previous section, we are assured that a (strictly positive) func- 
tion exists such that the vector field T = JT is complete. We look now for 
a reparametrization such that the vector field T coincides with the harmonic 
oscillator on the submanifolds of fixed negative energy. For this purpose, we 
notice that one can write T in terms of the energy in the following way: 

d / 1 u'V \ d 

T = Ua W^{^ ya - 2 ^ U lB^ (5 ' 14) 
so that, defining := {(y a , u a ) h £^ = E = cost}, one has: 

d f 1 n u PyP „\ d 
-En a - 2 — 



rL =u°— + 7 Ey a -2 — %^u a - — (5.15) 

dy a \2R 4 y R 2 J du a y ' 

This form of the field will soon prove useful. 

We have already noticed that, multiplying T, a second order vector field on 
M. = TM^, by a generic function / e ^(TM*), the reparametrized vector field 
is not second order on TR^; however, in the approach of par. 13.31 it is possible 
to introduce on M a new structure of tangent bundle with respect to which 
r is second order. One can consider as the starting functionally independent 
functions Y a = y a , so one has: 

C t (Y a )^fC r (y a )^fu a (5.16) 

Defining U a := f u a , one can build the following 1-1 tensor and vector field (as 
in Eq. I3~77jl : 

d d 

S:=dY a ®— — ; A ■=U a —— (5.17) 

dU a dU a 

so that: 

S(t) = A (5.18) 

that is, f 1 is second order. Its "second components" with respect to this new 
structure are: 

C t U a = fC T (fu a ) = fC T u a + fu a C T f (5.19) 

These general results make possible to reformulate our search for an / such that 
r coincides with an harmonic oscillator in the (simpler) search for an / such 
that F a are proportional to Y a , with a term possibly depending on the energy; 
more briefly: 

F a = g{£,L)Y a (5.20) 
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For simplicity we look first for a function only of the y a , that is / £ .F(R 4 ); so 
one gets: 

p a = p Fa + fuau (3jL f = 

If one chooses 4 / = 2R 2 , the two last terms cancel each other, and the first, 
when the energy is constant and negative, becomes the "second components" of 
an harmonic oscillator: 

f = 2R 2 -> F a = 2E C y a - 8i? 4 -^-w« + 8R 2 u /3 y l3 u a = ( 5 22 ) 
= 2£ £ y Q ' 

The reparametrized vector fields, in the new coordinates, reads: 

f = ^^ + 2 ( 2EV "i) yQ ^ (5 - 23) 
so, restricted to submanifolds of fixed energy = £ / T 1 (_E), it takes the form: 

d c) 
fL =U a - + 2EY a - (5.24) 

It is evident that this vector field, restricted to submanifolds of negative energy 
E, coincides with the vector field of an harmonic oscillator of frequency \J —2E: 
this was the relationship we expected. Put in different words, for negative 
energies we got a one-parameter family of harmonic oscillators, such that the 
frequency of each oscillator, and so its period, depend on this parameter, which 
is the energy of the conformal Kepler problem in R 4 ; when restricted to 
S , £-c is the pull back of the energy £,^ K of the Kepler problem in (see 
comment after Eq. In this way, the orbits lying on S with different 

have different periods, depending on £^ and so on £r_, K , thus can be related 
with the orbits of the Kepler problem, a non-isocronous system; this is what we 
expected taking into account the energy-period theorem. 

One can notice that the reparametrization involves the whole vector field, that 
is the eauation l5.23l is valid indipendently of the sign of £^ : so we have obtained 
that for — E > the reparametrized vector field coincides with a "repulsive" 
harmonic oscillator (that has open orbits, in agreement with the fact that the 
Kepler system has open orbits for positive energies), and for £^ = it coincides 
with the free particle. 

As we anticipated, in this way we have obtained the completion of the vector 
field T: it is possible to extend the submanifold S_e of fixed energy to include 
the points corresponding to the origin in the configuration space, the vector field 

4 the 2 is for convenience in calculations, it could have been any positive real number. 
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r being defined in these points since it coincides with the harmonic oscillator 
(attractive or repulsive) on each Y,e- In this way one gets a vector field, that we 
still call f, defined on the whole TR 3 , and here complete (cfr |2Z|, [2SJ)- 
The function / — 2R 2 , by wich one obtains the completion, is the very one 
whose existence the theorem we mentioned assures. 

An important comment is due about the kind of reparametrization we have 
choosen. We have pointed out that the function R 2 is one possible function 
that makes possible to build the relation between T and a family of harmonic 
oscillators, that is such that: 

F\^=g(E)Y a (5.25) 

it is by no means unique. One could have choosen, for example, R 2 multiplied 
by any other (positive) function of the energy, and so obtained as well: 

f= 2g(Z z )R 2 -> F a = 2g(Z L )Z L y a -8g(Z z )R 4 ^u a + &g(Zc)R 2 u (3 yPu a = 

= 2g(^)e L y a 

(5.26) 

that is on each £ e for E < a harmonic oscillator with frequency \f— 2 g(E) E. 
So, the choice of / = g(Hti)2R 2 brings just an overall multiplying factor in the 
expression of the reparametrized forces (reminding eauation l5.19|l : clearly the 
reparametrized vector field coincides, also in this case, with an harmonic oscil- 
lator (attractive or repulsive) on each S^, what changes is only the frequency 
of the oscillator (as in the case explicitely shown above). Our choice of / = 2 R 2 
(a very immediate one) was sufficient to achieve our purpose of relating T, and 
so the Kepler vector field, to a family of harmonic oscillators; one could use 
this arbitrariness in the choice of a multiplying function of the energy to require 
some additional properties to the reparametrized vector field. 
Moreover, we notice that the multiplication of a vector field in TR* by R 2 
correspond to the multiplication of his projection in TR 3 with respect to the 
Kustaanheimo-Stiefel map by r (see Eq. 14 .7|) ; actually, the reparametrization of 
the Kepler dynamical fields by the function f — r was a common tool to achieve 
its regularization (see e.g. |19) : what they actually did was the "changing 
of the time" from t to r such that dt = rdr). Gyorgyi (in ^3]) seems to be 
the only one who choose a different function / = r (in our notations); &c K 
corresponds to £^ (since the latter function is its pull back when restricted to 
S , the invariant submanifold of our reduction procedure), so this choice is in 
agreement with the ambiguity we have pointed out above. 

Summarizing, we have constructed an unfolding system T for the Kepler prob- 
lem and shown that it coincides, up to reparametrization, with a one-parameter 
family of harmonic oscillators whose period depends on the energy of the Kepler 
problem, such that when £^ < they are attractive oscillators (as we antici- 
pated on the basis of general considerations), when £^ > they are repulsive 
oscillators and the case £^ = corresponds to the free particle. For clarity we 
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briefly recall how the Kepler system can be obtained as a reduction of the un- 
folding system I\ One first restricts to the submanifold defined by the condition 
R 4 03 — (as we anticipated at the beginning of this section), that is invariant 
for r since Cr(R 4 0s) = 0. On this submanifold the action of the group U(l) 
given by the tangent of the right multiplication by exp(iA<73) provides an equiv- 
alence relation which is compatible with T. In this way one obtains the Kepler 
system as the reduced dynamical system. 

6 Constants of the motion 

In the preceding sections we have constructed an unfolding system for the Kepler 
problem and shown that, on submanifolds of fixed energy, it coincides, up to 
reparametrization, with a harmonic oscillator (attractive or repulsive). Using 
this last property, we can obtain the constants of the motion of this system 
starting with those of the harmonic oscillator. Then, we will characterize the 
one parameter group we use to perform the reduction, point out its Hamiltonian 
character and then obtain the subalgcbra of symmetry for the conformal Kepler 
problem in 4 dimensions that reduces to a symmetry algebra for the Kepler 
problem in 3 dimensions. 
For this purpose, we recall a known result. 

6.1 A preliminar result 

Let G be the Hamiltonian symmetry group of a given dynamical system, and q 
the corresponding symmetry algebra; let us consider a one parameter Hamilto- 
nian subgroup of G, whose Hamiltonian function is F and whose infinitesimal 
generator is Y. Let us reduce the dynamical system choosing the invariant sub- 
manifold S a = {m e M h F = a = const} and the equivalence relation on S a 
given by the one parameter group (which acts preserving E ): two points are 
equivalent if they are on the same orbit of the group. With these premises, if 
& is a subalgebra of g, we want to find under what conditions & is a symmetry 
algebra for the system obtained as reduction of the former by the action of this 
one-parameter subgroup. 

First of all, the functions fi <G g, restricted to S a , have to be constant on the 
orbits of the one parameter group, i.e. on the integral curves of the vector field 
Y: 

£ Y / = (6.1) 

(there is no need to explicitly restrict to E since the integral curves of Y belong 
to E ); since Y is Hamiltonian, this condition can be expressed by: 

{fi,F} = (6.2) 

This condition assures that the fi restricted to each submanifold E OJ i.e. on all 
the leaves of the foliation, can be obtained as the pull-back of a function on the 
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quotient manifold E a ; if one is interested in only one leaf, say E a , it is sufficient 
to require: 

{fi,F}h*=Q i-e- {fi,F} = c(F-a) (6.3) 
where c is a constant. 

Let us call fa the function on the reduced space E a that correspond to the fi, 
that is / o 7T = /. The condition 16. II for 16. 2|) assures that the /, are constants 
of the motion for the reduced system. Indeed, denoting $p the flow of T on E a 
and <f>|, the flow of T on E, it follows from the fact that T is the projection of 
r with respect to it: 

$j,0 7T = 7TO$^ (6.4) 

So one has, Vm 6 E Q , m = n(m) 6 E a : 

/(4(m)) = /(4 o 7r(m)) = /> o $*,(m)) = 

/o$*,(m) = /(m) = /(m) (6.5) 

i.e. / is a constant of the motion for the reduced system. 

Thus, the algebra of the constants of the motion for the reduced system is the 
subalgebra of the constants of the motion for the starting system that Poisson- 
commutes with the Hamiltonian for the one parameter group we use to perform 
the reduction. 



Pursuing the program we have stated, we will first characterize the symmetry 
of the unfolding system T and the one parameter subgroup we use to perform 
the reduction; we will notice that we are in the hypotesis of the previous result, 
which we will use to get the constants of the motion for the Kepler problem in 
3 dimensions. 



6.2 Digression: on the symmetry of the harmonic oscilla- 
tor 

In order to find the constants of the motion for the unfolding system T, we will 
study the ones of the reparametrized vector field T, since they are not affected 
by reparametrization. At this aim, our main tool will be its correspondence, on 
submanifolds of fixed energy, with an harmonic oscillator, since it allows us to 
obtain the symmetry group and the constants of the motion for T from the ones 
of the harmonic oscillator. 



We consider now an harmonic oscillator with four degrees of freedom, defined 
on TM 4 . We recall that it admits a Lagrangian description, with Lagrangian 
function Lho = 1/2 (U a U a — n 2 Y a Y a ), and corresponding Cartan 2-form: 

ojho = dY a A dU a (6.6) 

and related Poisson bracket, as in Eq. Ij3.8|l : 

if x = _dj__dg_ _ _dj__dg_ , 
\Ji9sho djj a dY a dY a dU a v D -'J 
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Introducing complex coordinates: 

z a = U a + inY a (6.8) 

the energy function reads: 

I 3 

the Cartan 2-form is: 

uho — (2m)~ 1 dz a A dz Q 
and gives the following Poisson Bracket: 

{f,9}HO = -2Kt(d a f-d a g-d a f-d a g); f,g e F(TR 4 ) (6.11) 

It follows that the quadratic functions that Poisson-commute with the Hamil- 
tonian arc of the form: 

F = (2Ki)- 1 c aP z a z (3 (6.12) 

where the c Q( 3 are costant complex numbers such that c a p — —c a /3, if we re- 
quire that constants of the motion are real. The matrix we associate with this 
quadratic form is a complex 4x4 antihermitian matrix C. Recalling that a 
complex antihermitian matrix can be decomposed in the form C = A + iB: 
A = —A T ,B = B T , A a p,B a p <G R with A and B real matrices, respectively 
antisymmetric and symmetric, one gets: 

F = A a fjL a p + BafjQaf) 

where L a p = (2n)^ 1 lmz a z^ 

Q afj = {2k)- 1 Rcz°'z /3 (6.13) 

with L a p — —Lp a and Q a p — Qp a ] moreover, if F = {2ni)^ 1 c a pz a z 13 , G = {2ni)~ 1 d a pz a z 13 : 

{F,G} = (2m)- 1 [C,DUz a z^ (6.14) 

So we have exstablished a correspondence (a Lie algebra isomorphism) between 
the quadratic constants of the motion and the complex antihermitian 4x4 
matrices: they form the Lie algebra u(4). 

Since one of the constants of the motion is the energy (the one corresponding to 
C = ikI), and since it is a central element, one usually identifies the Lie algebra 
su(4) as the symmetry algebra for the harmonic oscillator. The corresponding 
symmetry group is SU{A). 

6.3 Constants of the motion of the unfolding system 

Since T, restricted to submanifold of fixed negative coincides with a har- 
monic oscillator with frequency k = \/—2E, the constants of the motion for the 
harmonic oscillator are constants of the motions for T , restricting to XL , i.e. 



(6.9) 



(6.10) 
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the part of TRg with negative £^ ; since n is constant only on the submanifolds 
Te, hi the extension of the costants of the motion from each Te to £_ we have 
to replace it with \J — 

One obtains the following expressions 5 in the "new" coordinates (Y a , U a ): 
L at3 =\ {Y a U(3 - U a Y p ) 

Q«(3=\(U a Up-(2Z c )Y a Yp) (6.15) 

From this expression one can notice that L a p,Q a p are well defined indepen- 
dently of the sign of £,c , so can be extended to the whole space, not only to the 
portion where £,c < 0. Indeed, one can check that they are still constants of 
the motion for T, whatever the sign of £^ is. 

We point out that we haven't yet a Poisson structure for the constants of the 
motion of T, as we got for T. It seems a natural choice to use the one obtained 
from that of the harmonic oscillator. Since both the symplectic 2-form and 
the Poisson bracket of Eqs. i|6.6[l . (|6.7|l do not depend on the frequency of the 
harmonic oscillator, they can be extended naturally from each submanifold Tie, 
where T coincides with the oscillator of frequency \/—2E, to the whole space: 

= dY a A dU a 
r, _dj^_dg_ _ _dj^_dg_ 

U,9t~ dU a dY a dY a dU a { ■ ' 

In this way, the space of the constants of the motions for T becomes a Lie algebra 
under the Poisson bracket 16.71 which we denote Tf., and to each function / in 
Tf we may associate a the vector field X/ such that ix/W^ = d/. 
As we have already pointed out, the constants of the motion do not depend on 
the parametrization, so L a p and Q a p are constants of the motion also for T. 
We notice that, expressing in the "old" coordinates (y a , u a ), comparing with 
Eq. 15.121 one has: 

1 

w~ = -jUz, (6-17) 

i.e. they are the same, apart from a constant factor that can be absorbed. This 
makes possible to build a Lie algebra isomorphism between the algebras J-f of 
the constants of the motion for T and J-f of those of T. 

As for the expression of the Poisson bracket we got in terms of z and z, if we 
replace k with \/—2E in it, we would violate the Jacobi identity (note that 
Eqs. I|6. 8JI do not define a changing of coordinates when k is replaced with 
\/—2E). If, however, we restrict ourselves to functions that Poisson commute 
with £^, we would again satisfy the Jacoby identity on this restricted subalge- 
bra, because on it En is a central element. This remark will prove useful in a 
following section, where we will be allowed to use Eqs. (|6.11|) and l|6.14|l since 
we will deal only with constants of the motion. 

5 we have rescaled the Q a p by a factor y / —2£,£ ; it is a central element of the algebra of the 
constants of the motion we will introduce later 
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With these premises, we are able to investigate in a better way the one pa- 
rameter group we perform the reduction by and then get its expression in the 
"new" coordinates (Y a , U a ) and find the associated constant of the motion. As 
we have already stated, this group acts as the tangent lift (with respect to the 
"old" tangent bundle structure on TR*) of the right multiplication by exp(iAcr 3 ): 



gexp(i\cr 3 ) 



or, in coordinates (y a ,u a ): 



( y 1 \ 



V y° J 



( cos A — sin A 

sin A cos A 

cos A — sin A 

y sin A cos A J 



( y 1 \ 

v 2 
y 3 

\y° J 



(6.18) 



(6.19) 



or, calling S\ the 4x4 matrix above and y* the column vector in . 



(6.20) 



Its infinitesimal generator is the left invariant vector field usually denoted as 
X 3 : 

„ d , d , d „ d 

X3 = y i^-y 7^ +y i^-y ^ 



dy 3 » dy° * dy 2 
The tangent lift of the above action is given by: 

(y*,u*)^(S A y*,S A u*) 

and its infinitesimal generator is the tangent lift of X 3 : 



(6.21) 



(6.22) 



d 



d 



d 



d 



x 3 =^v _2/ V +y V _2/ V 



d , d , d 2 d 

du 3 du° du 2 du 1 



This vector field is Hamiltonian with respect to the 2-form u>c and its Hamil- 
tonian function is given by: 

since X^ is the tangent lift of X 3 and ojc — — d#£, one obtains: 



h = AR% 
which in (y a ,u a ) coordinates reads: 

h = 4R 2 (y°u 3 - y 3 u° + y V - yV) 



(6.24) 
(6.25) 



(6.26) 
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If h is a constant of the motion for T (as it should be, in order to perform the 
reduction procedure we stated), it is a constant of the motion also for T. 
In the "natural" coordinates for the new tangent bundle structure one has: 



h = 2{Y°U 3 - Y 3 U a + Y X U 2 
Comparing with Eq. (|6.15(1 . we notice that: 

h = 4(L 30 + L 12 ) 



Y 2 U i ] 



(6.27) 



(6.28) 



i.e. h is the constant of the motion for T to which corresponds the antisymmetric 
3x3 real matrix: 



N 3 



( 


-i 








\ 


1 




















-1 




V o 





1 








(6.29) 



The vector field associated to this constant of the motion by means of the 
symplectic structure lo^ is: 



X = Y { 



8 

+u° 



d 



Y 1 



dU 3 



dU° 



d 

u 1 



Y< 



dU 2 



d 

W 1 



dU 1 



(6.30) 



Summarizing, we are in the position to use the result of section l6~T1 so, to char- 
acterize the symmetry algebra of the reduced system, we only have to find the 
constants of the motion which satisfy Eq. H6.3(l with the function h Hamiltonian 
for X3, i.e. {fi, h} — c ■ h (since ft = on So), where c=const. 
For simplicity we first consider the case c = 0, and then show that the functions 
obtained in this way are all the constants of the motion for the reduced system. 
For this purpose we can use the isomorphism between the Lie algebras of the 
quadratic constants of the motion and of the complex 4x4 antihermitian trace- 
less matrices, previously extablished: according to Eq. (|6.14(l . a quadratic func- 
tion that Poisson-commutes with h corresponds to a matrix that commutes 
with N3. Recalling that C = A + iB; A = —A T , B = B T , one can write 
the most general real traceless symmetric (respectively antisymmetric) matrix 
and then impose that its commutator with N3 vanishes. After a lenghty but 
straightforward calculation one obtains the following basis for the Lie algebra 
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that commutes with N3 (cfr. |3(J|1: 



2Mi = 



2AU = 



2D 2 = 



( 



V 1 



1 \ 



V 



1 



2M 2 = 



2L>i = 



2D* = 



i 

i 

V 



(6.31) 



where the Mj are real antisymmetric matrices and the Dj are symmetric pure 
imaginary ones obtained as Dj = iBj, with Bj real symmetric. One can check 
the following commutation relations among these matrices: 

[Mi, Mj] = e ijk M k 
[Di,Dj] = e ijk M k 

[M i ,D j }=e i j k D k (6.32) 

These elements close a subalgebra of su(4); defining A{ = h(Mi + Di) and 
Bj = \[Mi — Di), from their commutation realtions: 

[Ai, Aj] — €ij k A k 
[Bi, Bj] = €ij k B k 

[Ai,Bj]=0 (6.33) 
one recognizes the algebra su(2) x su(2). 

So we have obtained the result that the symmetry subalgebra that commutes 
with h is isomorphic to su(2) x su(2). 

We can now express the constants of the motion corresponding to the Mj , Dj 
as functions on TM 4 using Eqs. (|6~T2j) . lj(H5|l : 



J l = 


1 

2 


(r 1 ^ 




V Y 3 U 2 


- Y 2 U 3 


J2 = 


1 

2 




- Y 3 U X - 


\-Y 2 U° 


~Y a U 2 


■h = 


1 

2 


{Y X U 2 


- Y 2 U X - 


h y0 [/ 3 


-Y 3 U° 



Qx = - {V X V 3 + U 2 U° - 2E L {Y 1 Y 3 + Y 2 Y )) 
Q'2 = ^ (U 2 U 3 - U 1 ^ - 2£ C (Y 2 Y 3 - y 1 ^ )) 



(6.34) 



h = ~ [(U 1 ) 2 + (U 2 ) 2 - (U 3 ) 2 - (U ) 2 - 2^((Y 1 ) 2 + (Y 2 ) 2 - (Y 3 ) 2 - (Y ) 2 )] 
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We notice that, since the last three constants of the motion depend on the 
energy £r_, when we extend them from the submanifold Eg to the whole space, 
the factor y/— 2£^ is no more constant: the Poisson commutation relations 
become: 

{Jit J j} ^ — ^ijkJk 
\JiiQj\ r*j — ^ijkQk 

{QuQj}~=e ijk Jk(-2£z) (6.35) 

Since we are interested in the subalgebra of the constants of the motions for T 
that correspond to the symmetry algebra for the Kepler problem, we now express 
the Ji and Qi in the coordinates (y a , u a ), since they are the natural coordinates 
for the structure of tangent bundle with respect to which T is second order: 

- y 3 u 2 — y 2 u 3 ) 

- y 2 u° - y°u 2 ) 
- y°u 3 - y 3 u°) 

ZcivV-yV) (6.36) 
i ££ (( 2 / 1 ) 2 + (2/ 2 ) 2 -(y 3 ) 2 "(y°) 2 ) 

Recalling Eq. (|6.17|) , we have to rescale these functions to get the same commu- 
tation relations of Eq. (|6.3[) with the Poisson bracket { , }x ■ 
After this remark, we can recognize in the Ji,Qi the functions in .F(TRq) cor- 
responding to the angular momentum and the Runge-Lenz vector in JF(TM[|) 
respectively. Since is a central element, also in this case we can rescale the 
Runge-Lenz vector, so to obtain an so(4) algebra when £^ < and an su(2) 
when £ £ > (cfr. Eq. lPHU|> and below). 

This last consideration was sufficient to conclude that the functions that Poisson 
commute with h are all the constants of the motion for the reduced system; so 
it is not necessary to pursue the investigation also for values of c different from 
0. One can arrive at the same conclusion noticing that the above functions (ac- 
tually five of them) are independent and generate the whole space of constants 
of the motion for the Kepler problem, its whole symmetry algebra. 
The constants of the motion we have obtained coincide with those of how- 
ever we feel that our Lagrangian derivation clarifies the way they arise without 
using any educated guess. 



1 fi It"- 



Ji = ^- 2 (y 1 u°-y°uU 
■h = ^{y 1 u 3 -y 3 u 1 ^ 

J 3 = i ^(y 1 « 2 -yV4 
Qi = g^tuV + uV)- 
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7 Conclusions 



The leading idea of this paper was that completely integrable systems should 
arise from reduction of linear ones. We have considered as case of study the 
Kepler problem in three dimensions. After consideration based on general prop- 
erties of the system, we have arrived at a possible unfolding system, recovering 
the known relation of the Kepler problem with a family of harmonic oscillators, 
in a constuctive way, as far as it has been possible. Moreover, we have charac- 
terized the symmetry algebra of the unfolding system and the subalgebra of the 
constants of the motion for the Kepler system. 

It seems to us that the procedure we presented could be useful in the unfolding 
of general supcrintegrable systems. 

Following the same line we hope to develop an unfolding procedure also for 
quantum systems. 
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